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DIFFERENCE EQUATIONS AND THEIR APPLICATIONS 


Louis A. Pipes 


I. Introduction. 

During the last fifteen years there has been a great deal of activity 
in the applications of difference equations in the solution of many prob- 
lems that arise in statistics, science, and engineering. 

The development of high-speed digital computing machinery has moti- 
vated the use of difference equations as approximations to ordinary and 
partial differential equations. [See bibliography, 1, 2, 3]. The use of dif- 
ference equations for solving partial differential equations was discussed 
in 1928 in a celebrated paper by Courant, Friedrichs, and Lewy. [4]. How- 
ever, the methods presented in this paper were not put into practical use 
until about 1943 as a result of the stimulus of war-time technology and 
with the assistance of the first digital automatic computers. Problems in- 
volving time-dependent fluid flows, neutron diffusion and transport, ra- 
diation flow, thermo-neuclar reactions, and problems involving the solu- 
tion of several simultaneous partial differential equations are being solved 
by the use of difference equations throughout the country. 

Besides the use of difference equations as approximations to ordi- 
nary and partial differential equations, they afford a powerful method for 
the analysis of electrical, mechanical, thermal, and other systems in 
which there is a recurrence of identical sections. The study of the be- 
havior of electric-wave filters, multistage amplifiers, magnetic ampli- 
fiers, insulator strings, continuous beams of equal span, crankshafts of 
multicylinder engines, acoustical filters, etc., is greatly facilitated by 
the use of difference equations. The usual methods for solving such sys- 
tems are generally very lengthy when the number of elements involved is 
large. The use of difference equations greatly reduces the complexity 
and labor in problems of this type. 

II. Nature of the Calculus of Finite Differences. 

As is well-known, the most important concept of mathematical an- 
alysis is that of a function. If to a given value of x a certain value of y 
corresponds, we say that y is a function of the independent variable r 
and write symbolically y(z): Two types of functions occur in the appli- 
cations of mathematical analysis to physical problems; first there are 
functions in which the variable x may take every possible value in a given 
interval, that is, the variable z is continuous. Functions of this type are 
studied in the ordinary Differential and Integral Calculus. 
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A second type of function exists in which the variable z takes only 
the given values z,, 2, %,, -» 2,. In this case the variable z is discon- 
tinuous. The methods of the ordinary calculus are not in general appli- 
cable to such functions. The Calculus of Finite Differences deals es- 
pecially with functions of the second type, but it may be applied to func- 
tions of the first type as well. The origin of the Calculus of Finite Dif- 
ferences is usually ascribed to Brooks Taylor who published his, “Me- 
thodus Incrementorum” in 1717. He pointed out that the relations between 
the terms of arithmetic and geometric progressions, which were known for 
a long time represented the simplest examples of difference equations. 

There are many analogies between the Differential and Integral Cal- 
culus and the Calculus of Finite Differences. For example, in the Cal- 
culus of Finite Differences, the first difference, denoted by A(A)F(z) is 


defined by, 
(2.1) A(A)F (2) = F(x+h) - F(a) 


In the Differential Calculus, the first derivative of a function F(x) 
is defined by the equation, 


(2.2) DF(z) = Limit _ Limit 


h 


It can be shown that the ath. derivative of /'(z) may be expressed in 
the form, 


(2.3) = Limit 


The operator A(A) is called the forward difference operator. When 
there cannot be any misunderstanding about the increment A, the differ- 
ence operator A(h) is merely written as A, 


Basic Operators used in the Calculus of Finite Differences. 


The following fundamental operators are used widely in the Calculus 
of Finite Differences. 
1. The Displacement Operator E. 


(2.4) EF(2) = F(a+h) 

2. The Forward Difference Operator A. 

(2.5) AF (a) = F(x+h) -F (a) 
3. The Backward Difference Operator V . 

(2.6) VF (2) = F(z) -F(a2-h) 
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4. The Central Difference Operator 5. 
5F (x) = F(a+5) F(z 5) 
5. The Averaging Operator yu. 
(2.8) (a) = + 
6. The Symbolic Representation of Taylor’s Expansion. 


(2.9) F(z) = Fla+h), 
dz 


Equivalent Operators. 
By using the above fundamental definitions of the basic operators 


used in the calculus of finite differences, the following equivalent opera- 
tors may be easily deduced. 


(2.10) 
(2.11) A = (E-1) 
(2.12) V =(1-E7') 
(2.13) 


(2.14) 


It can be shown [5, 6, 7, 8] that except for a few restrictions, these 
operators can be manipulated according to the laws of ordinary algebra, 
and many relations between the basic operators may be obtained. 

For example, we may write, 


(2.15) | D= tale by the use of (2.10) 


(2.16) 
(2.17) 5 = 2sinh(hD/2) 
(2.18) pw = cosh(hD /2) 


An interesting and useful exposition of the symbolic calculus of op- 
erators will be found in a paper by Bickley [9]. By the use of these oper- 
ators, and their relations between them, several finite difference formulae 
for interpolation, differentiation, integration, scanning, etc. may be ob- 
tained. These formulae are of great importance in numerical analysis. 

For example, if the logarithm in (2.15) is expanded in powers of A, 
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we obtain, 


(2.19) p = _ AP_A*, 
h h 2 3 4 

This expression may be used to compute the derivative of a function 
in terms of its forward differences. 
The Scanning Operator. S(h) 

In the theory of scanning, if we are given a function F(z), it is re- 
quired to obtain the function (2) defined by, 

2) 

(2.20) = | Fladdz 


(t—h; 2) 


If (2.20) is differentiated with respect to z, the result is, 


Dela) = 
2 


= 5F (2) 


= F(a) = S(A)F (2) 


Where S(A) is the scanning operator defined by, 


7.99) SA) 5 2sinh(hD/2) hsinh(AD /2) 
D D (AD /2) 


If d(x) is known and it is desired to determine F(x), we may write, 
2.24) F(x) = S~ A) A(z) 


In (2.24), S~*(h) is the inverse scanning operator. 


III. Linear Difference Equations with Constant Coefficients. 

The most important difference equations that arise in applications 
to dynamics, electrical circuits, structural problems, etc., are linear dif- 
ference equations with constant coefficients. 

These equations are of the type, 


In the usual applications of difference equations, h = 1, so that the 
operator F is defined by, 


E*y(2) = ylat+k) 


Hence, 
(3.2) 
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Difference equations are classified in the same manner as differen- 
tial equations. The equation (3.1) is called a linear difference equation 


of order n with constant coefficients, provided that the coefficients a, 
are constants. If V(x) = 0, then the equation is said to be homogeneous; 
if not, it is said to be non-homogeneous. 

Let, 


(3.3) P(E) = (a,E"+a,_ 


The general solution of (3.1) may be written in the form, 


(3.4) y = + 


If this is substituted into (3.1), the result is, 
(3.5) P(E) =0 


Let, 
(3.6) P(E)y (x) = 0, 


and 


(3.7) P(E = V(z) 


By analogy to the solution of linear differential equations with con- 
stant coefficients, y, is called the complementary function, and Vp the 


particular integral of the differential equation (3.1). 
Determination .of the Complementary Function. 

To determine the complementary function y,(z) we assume a tenta- 
tive solution of the form, 


(3.8) y = Aq” 


where A is an arbitrary constant and g a number to be determined so that 
(3.8) is a solution of (3.6). It is noticed that if the operator EF is applied 
to (3.8) the result is, 


(3.9) (2) = E*Ag* Aq? q*y (a), 


It is therefore evident that if (3.8) is substituted into (3.6), the re- 
sult may be written in the form, 


(3.10) P(E)y (x) = P(g)y,(x) = 0 


If the trivial solution y (2) = 0 is excluded, then we must have, 
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(3.11) P(g) = 0 


This is the characteristic equation of the difference equation (3.1). 
The characteristic equation is an algebraic equation of the nth. degree 
in g. The roots of this equation determine the possible values of g that 
allow (3.8) to be a solution of (3.6). 

If the roots of (3.11) are the distinct numbers q,,-+-q,, then the gen- 
eral solution of (3.6), may be written in the form, 


(3.12) = 


where A,, A,,--- A, are arbitrary constants. 

This solution is valid provided the roots g; are distinct, regardless 
of whether they are real or complex. 

If the characteristic equation (3.11) has a multiple root q, of the rth. 
order, it can be shown [7] that the part of the complementary function y_, 


that corresponds to this repeated root is, ' 
(3.13) Yo 2) = (C Cyr? 


where the C;’s are arbitrary constants. 
The Particular Integral, y 

The differential equations most frequently encountered in practice 
are linear and homogeneous with constant coefficients. The various meth- 
ods for obtaining the particular integrals of linear equations with constant 
coefficients have their counterpart in the theory of difference equations. 
[5, 6, 7, 8]. 

The method to be used in determining the particular integral Vp de- 
pends on the nature of the function V(x). To illustrate the general proce- 
dure, let the function V(x) be of the form, 


(3.14) Vir) = Ke™* 
The particular integral Vy in this case satisfies, 
(3.15) P(E)y, = Ke™* 


In this case, to obtain the particular integral assume a solution of 
the form, 


(3.16) 


where 2 is to be determined. It is noticed that, 


(3.17) E*Be™® Ae™* emt 
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It is therefore evident that, 
(3.18) P(E)Be™® = Ae™*P(e™) 
Hence if (3.16) is substituted into (3.15) the result is, 
(3.19) Be™*P(e™) = Ke™* 
Therefore if P(e”) 4 0, B is given by, 
(3.20) B=K/P(e™), 


and the particular integral of (3.15) is 


(3.21) yp = Ke™*/P(e™) 


In general, the particular integral of (3.1) may be expressed symboli- 
cally in the from, 

(3.22) = P~"(E)V (2) 

If V(x) is a rational and integral function of z it is possible to ex-, 
pand P~‘(E£) in a series of ascending powers of A by the-substitution 
E = (1+) and operate with this series on V(z). If V(z) is a rational and 
integral function of x, the result of this operation will give a finite num- 
ber of terms. 

In some cases it is convenient to expand the operator P~‘(E) into 
partial fractions. If this is done, the required particular integral may be 
obtained in terms of indefinite sums. [7]. 


IV. Applications of Difference Equations to the Analysis of a Uniform 
Electric Network. 

An interesting and typical technical problem that illustrates the use 
of difference equations, is to determine the current and potential distri- 
bution of the electric circuit shown in Figure 1, in the periodic steady 
alternating state. 


— 


Figure 1. 
The circuit of Figure 1 is a cascade connection of n boxes. Each 
box contains the Pi Network of impedances Z and admittances Y/2, shown. 
in Figure 2. 
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Circuit Inside the kth. Box 
Figure 2. 
Each box of the network of Figure 1 contains the Pi Network of Figure 2. 
This consists of a symmetrical arrangement of two shunt admittances 
Y/2 and a series impedance Z. 

If the sending end potential V, and the sending end current /, are 
specified, it is required to determine the potential and current distribution 
of the uniform network. If Kirchhoff’s first and second laws are applied 
to the network of Figure 2, the following simultaneous difference equa- 
tions are obtained: 


1 


In order to solve the simultaneous system of difference equations 
(4.1) and (4.2), it is necessary to separate the variables. If (4.2) is solved 


for /,_,, the result is, 
1 Y 
If this equation is used to eliminate /,_ , and /, in (4.1), the result is, 


(4.4) Ving = 0 


This is a homogeneous difference equation of the second order of the 
type discussed in Section III. Its solution is of the form, 


(4.5): V, = 


where A is an arbitrary constant and g satisfies the characteristic equa- 
tion, 


(4.6) q?-(2+ZY)q+1 =0 


The two roots of this equation, satisfy the equations, 
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=1 


= (2+ZY) 


then 

(4.10) 91+ = (2+ZY) = 2cosh(a) 
The Propagation Constant of the Network. 


(4.11) a = cosh~ '(1+ZY/2) = Propagation Constant. 


(4.12) Vi 


Z 


The Characteristic Impedance of the Network. = Z, 


(4.14) Z, = tanh(a/2) 


2 
Y 
Determination of the Arbitrary Constants. 


The two arbitrary constants 4, and A, may be obtained from the 
boundary conditions. 


(4.15) V,=A,+A, 


l 
( 4.16) = —(-A,+A,) 
“0 
Hence, 


(4.17) A, = 
(4.18) A, 


If these values of A, and A, are substituted into (1.12) and (4.13) 
the result may be expressed in the form, 


(4.19) V, = V ,cosh(ka) Z sinh(ka) 


(4.20) I,=-V, ],cosh(ka) 
“0 
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cosh(a) Z ,sinh(a) 
(4.21) [T] = = Transmission Matrix. 


sinh(a) 


Z, 


Vi 
(4.21) (y);, = 
I, 


Then, 
(4.22) (y), = [TI-*(y), 


cosh(a) 


and 


(4.23) = TY), 


V. Longitudinal Oscillations of a Long Uniform Train. 

Another typical example of the use of difference equations is the 
application to the dynamical problem of determining the resonant fre- 
quencies and modes ofoscillations ofa chain of n cars coupled by springs 
that satisfy hooke’s law as shown in the Figure 3. 


| 


4 


Q OO 690 
? 


Idealized Uniform Train of n Cars 
Figure 3. 

Figure 3 represents an idealized train of n cars of equal mass m 
coupled to each other by identical springs of spring constant k. The ef- 
fects of friction are neglected. The simplest method of analyzing the be- 
havior of this system is by the electrical analog method which makes use 
of the results of Section IV. [See 10 in the Bibliography] 

In this section, the dynamical system of Figure 3 will be analyzed 
by the use of a chain of differential equations and difference equations. 

If Newton’s second law of motion is applied to the rth. car in the 
chain of cars of Figure 3, the following differential equation is obtained, 


(5.1) mar’ + ,)+ = 0, r= 2, 3, +++ (n—1) 


In order to study the oscillations of the system, assume a solution 


of (5.1) of the form, 
(5.2) = y,sin(wt+6) 
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where y, is an amplitude function and is an angular frequency to be 


‘ determined. 6 is an arbitrary constant. 


If (5.2) is substituted into (5.1) and the result divided by sin(wt+8), 
one obtains, 


(5.3) = 0 
where 
(5.4) 6 = (1-w?m/2k) 


To solve the difference equation (5.3) assume, 
(5.5) y, = Aq’ 
If this solution is substituted into (5.3) the following characteristic equa- 
tion is obtained : 


(5.6) q?-20g+1=0 
Let the two roots of this equation be g, and q, and let, 
(5.7) Jaze 


Therefore, 
(5.8) (9 ,+94)/2 = (e%+e~%)/2 = cosh(a) = 0 


The solution of the difference equation (5.3) may be written in the 
convenient form, 


(5.9) y, = B ,cosh(ar) +B,sinh(ar), r= 2,3, 4, ++ (n-1) 


where B and B, are arbitrary constants. 
The Boundary Conditions. 

In order to obtain the constants B, and B,, use must be made of the 
boundary conditions at the front and rear of the train. It is realized that 
there are no cars corresponding to r = 0 and to r = (n+1). In order to ex- 
tend the solution (5.9) to take into account r = 0 and r = (n+1) assume 
that there are phantom cars that correspond to these indices. Since the 
front and rear of the train is able to move freelii, there cannot be any 
strain in the phantom springs connecting the car r = 0 to r=1 or the spring 
that is supposed to connect the car r =n to the car r = (n+1). This condi- 
tion may be expressed mathematically by requiring that, 


(5.10) Yo=41 
(5.11) 


and (5.9) can be used as if it were valid for r = 0 and r = (n+1). 
The boundary conditions (5.10) and (5.11) lead to the two equations. 
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(5.12) B = B ,cosh(a) +B ,sinh(a) 
(5.13) B ,cosh(an)+B,sinh(an) = B ,cosh{a(n+1)]+B ,sinh[a(n+1)] 


B, can be expressed in terms of B, by (5.12) in the form, 
(5.14) B, = -tanh(a/2)B , 


If this expression for B, is substituted into (5.13) and some algebraic 
reductions are made, the result is, 


(5.15) B jtanh(a/2)sinh(an) =0 


For a non-trivial solution, we must have B, = 0. Hence the possible 
values of a must satisfy the transcendental equation, 


(5.16) tanh(a/2)sinh(an) = 0 


This equation is satisfied provided that, 


(5.17) sinh(an) = 0, 


or 
(5.18)  a@=a,= jas/n, 8 =0,1,2,3,+-(n-1). j= (-1)* 


If s is given the value s = n, (5.16) leads to an indeterminate form 
and it is not satisfied. 
The Resonant or Natural Frequencies of the Train. 

To determine the possible natura] angular frequencies of the system, 
we may use (5.8) written in the form, 


(5.19) cosh(a) = cosh(a,) = = (1-w?m/2k) = cosh(jzs/n) 


Hence, 
(5.20) 2m/2k = 1-cosh(jzs/n) = 1-cos(zs/n) = 2sin*(7s/2n) 


Therefore, 
(5.21) w, = 2(k/m)*sin(ms/2n), = 0, 1, 2, 3, «+ (n=) 


These are the n natural angular frequencies of the principal modes 
of oscillation of the train. The frequency w, = 0, corresponds to a rigid 
body translation of the entire train. 

The Natural Modes of Oscillation of the Train. 

To determine the natural modes of oscillation of the system, write 

(5.9) in the form, 


(5.22) y, = B ,cosh(ar) + B ,sinh(ar) 


= B [cosh(ar)-tanh(a/2)sinh(ar)], r= 1, 2,3, 
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To each value of a = a, = jms/n, there corresponds an amplitude func- 
tion, 


(5.23) = 8, 8), 
where 
(5.24) f(r, s) = 


r=1,2,3,-n (cars) 
s = 0,1,2,---(n-1) (modes) 


The general solution of the chain of differential equations (5.1) with 
the appropriate boundary conditions can now be obtained by substituting 
(5.23) into (5.2) and summing over all the possible modes s. 

The general solution of the chain of equations (5.1), is therefore, 


s=(n-1) 
(5.25) 2, = B ,t+8,) 
4 n 2n n 


The index r refers to the car number and the index s pertains to the 
mode. The quantities B, and 6, are the 2n arbitrary constants of the gen- 
eral solution. They may be obtained if the initial displacements z,(0) and 
the initial velocities z/(0) of each car is given. 


VI. Difference Equation Method for Solving Partial Differential and Ordi- 
nary Differential Equations. 

There exists a vast literature in which the solution of ordinary and 
differential equations may be solved by the device of transforming the 
differential equation into an approximate difference equation. (1, 2, 3, 11] 

From the definition of the derivative operator D and the operators A, 
V, and 6, it is easily seen that the following operators are approximations 
for the derivative operator D, provided A is sufficiently small. As A+0 the 
right members of these equations tend to D. 


(6.1) 
pz”. 
(6.2) 

5 
=. 
(6.3) 


' To illustrate the use of these approximations to the operator D, let 
it be required to solve the differential equation, 


2 
(6.4) =2, or D*y=z. 
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Let the initial conditions be 
(6.5) y(0)=0, and Dy,0)=1 
For a sufficiently small 4, we have approximately, 


2 -1 
Hence the differential equation (6.4) may be replaced by the approxi- 
mate difference equation, 


h? 


The initial conditions can be replaced by, 


(6.8) (0) = 0, 


Let the difference equation be required to be satisfied at the successive 
points, z, = A, x, = 2h, ---.2, = kh. Then the difference equation (6.7) 


may be written in the form, 
(6.9) View + (k= 1,2, +++) 
where 
(6.10) y, = yla,) = ylkh). 
The initial conditions take the form, 
(6.11) Yo =9 
The solution of (6.9) is 
(6.12) = (h-h®/6)k + 
This corresponds to, 


2 
(6.13) y(r,) 2, 


The exact solution is, 
(6.14) = 2+ 2°/6 


If h is so small that h? can be neglected, (6.13) tends to (6.14). 
Partial Differential Equations. 

The Poisson differential equation will be taken as an example of 
the use of difference equations to solve partial differential equations. 
[3]. This expation has the form, 


(6.15) (Di+D ule, y) = fla, y) 


In this case we have approximately, 
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(6.16) D2 = 


2 
(6.17) D, = )/h 


where, 
(6.18) E ula, y) = ula+h, y), E ule, y) = ula, y+h) 


If the operators (6.16) and (6.17) are used, the partial differential 
equation (6.15) may be written in the approximate form, 


(6.19) wlr+h,y) y)+ y) + ulz,y-h) = 


The solution of this equation with prescribed boundary conditions is 
solved by Rosenbloom [3]. 


VII. Additional Applications. 

Additional applications of Difference equations have been made to 
the study of intermittent control systems [12] and the analysis of systems 
with hereditary characteristics [13], as well as in certain problems in the 
theory of probability [14]. 
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Editor, MATHEMATICS MAGAZINE, — 

Success to Professor D. H. Hyers, and praise for his establishment of 
the S.I.M.N.T.! (See Mathematics Magazine, -March-April, 1959.) Please 
enroll me as a member, and for my Habilitftationvorschlag I propose to 
replace “denumerable” and “enumerable” (which sounds too much like 
“innumerable”) by “countable”. 

Let the Executive Committee also take under consideration words 
whose primary meanings do not suggest their mathematica! significance, 
such as “field”, “ring”, “loop”, “histogram”, “regression”, etc. 


And can we get rid of the word “imaginary”, so misleading to-day, 
and substitute the single letter, “:” (or “j”), both as noun and as adjec- 
tive, with the definition i? = -1, and drop the -1 ? 

William R. Ransom 
Tufts University 


EULERIAN NUMBERS AND POLYNOMIALS 


L. Carlitz 


1. Introduction. Following Euler (6, pp. 487-491], we may put 


(1.1) .. A#1), 
(A#1) 


where H, = H,[A] is a rational function of A; indeed 


(1.2) R= [A] = Q-D"H_ 
is a polynomial in A of degree n-1 with integral coefficients. If we put 


nr 


(1.3) R= 2. A, (n> 1), 


s=] 


then the first few values of A,. are given by the following table, where 
n denotes the row and s the column; 

1 

l 


1 


7 


Alternatively, Worpitzky [15] showed that the A, may be defined by 
means of 


n 
(1.5) 
s=] 


The rational functions H, were studied in great detail by Frobenius 

[7], who was particularly interested in their relationship to the Bernoulli 

numbers. More recently Vandiver [14] has also made use of this relation- 

ship to obtain new properties of the Bernoulli numbers. Other recent oc- 

currences are [1], [2], [12]; generalizations occur in [4], [5], [13]. In view 
247 


(1.4) 
ll ll 
26 ~=—s-: 66 26 1 
57 302 302 
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of the long history of H, and A, , it is rather curious that, on the whole, 
these quantities are not very well known. Indeed an examination of Math- 
ematical Reviews for the past ten years will indicate that they have been 
frequently rediscovered. Actually there is no detailed discussion of H, 
in any book. On the other hand, Riordan, in his recent book [11], does 


develop a few basic properties and indicates the connection of A, , with 


certain combinational problems. 
The present paper is mainly expository. We include numerous proper- 


ties of H,, and the related polynomial 


n 
(1.6) H (ula) = 

r=0 
indicate the connection with Bernoulli numbers and polynomials and final- 
ly obtain some arithmetic properties of H,. For the combinatorial appli- 
cations the reader is referred to Riordan’s book [11]. The #/,, also occur 
in certain criteria for Fermat’s last theorem; this is discussed at length 


in Bachmann’s book [3] and will not be considered in the present paper. 


2. The defining relation (1.1) is evidently equivalent to 
(2.1) (H+1)"=AH, (n>0), Hy=1, 


where after expansion of the left member, H” is replaced by H,; we shall 


use this convention frequently. If f(z) is an arbitrary polynomial in 2, 
(2.1) implies 


(2.2) A H+1) = Af(H) + 


In particular, for f(z) = (7), we get 


(2.3) (Ast)  (m20, 


which implies 
H H 
(2.4) (m2. 
Repeated application of (2.3) gives 
H H 
 (m20); 


in particular we have 


(2.5) ) 4. 
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It also follows from (2.2) and (2.3) that 
H H 
(2.6) = (m2); 


‘in particular by (2.5) 


(2.7 wig 


Again, if f(z) is an arbitrary polynomial of degree n, we recall that 


n 
fe) =) 


where 
Aflz) = fla+1)-flz), 2) = f(2+1)-A™' 
Using (2.5) we get ; 


n 
(2.8) fH) = > 
r=0 


Since 


r 


ATA0) = 


we may write 


n 


It is perhaps of interest to mention that (2.9) can also be obtained as 
follows from (1.1). For |z| sufficiently small we have 


(A-1)7"(e7-1)" 


r=0 


n 


s=0 
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r=0 
r=0 s=0 
= 
( 
: |__| 
n 
r=0 
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so that 


n 
r=0 s=0 


and (2.9) follows at once. 
Returning again to (1.1) it is easily verified that 


(2.11) = 
also differentiation yields 
d 
(2.12) (n+ DAR, + 
where FR, is defined by (1.2). From (2.11) and (2.12) we immediately ob- 
tain 
(2.13) A. = 
and 


(2.14) A = sA , + (n—8+2)A, .. 


n+1,8 n,s— 


By means of (2.14) one can easily extend the table (1.4). A convenient 
check is furnished by the formula 


n 
= > (n>1). 
s=1 


We also note that 2, [-1] = 0 for n>1. 
Frobenius remarks that it follows from (2.12) that the n-1 roots of 
RA] = 0 are real, negative and distinct; also for each root A, the recip- 


rocal A, is also a root. Moreover the roots of F,,, ,[A] = 0 are separated 


n+i1 
by the roots of R, [A] = 0. In the next place, by (2.10) 


n 
r=0 j=0 


n 


r 


r=0 s=0 j=0 
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n 


>: As 


s=0 7=0 


s=0 j=0 
Hence, by (1.3) and (2.13), we get 


8 
(2,15) Ang = >, (DAE 
j-0 


a formula due to Euler. 


n n 
r=0 r=0 


Since 


it follows from (2.10) and (2.15) 


n n 
s=] s=] 
which establishes (1.5). This proof is taken from Frobenius. 
3. We now consider the polynomial 


(3.1) = H,(u| A) = (u+H)" 


defined by (1.6). We evidently have the generating function 


(3.2) 


n=0 


It follows from (3.2) that 
(3.3) H,(u+1) -AH,(u) = (1-A)u” 


Moreover (3.3) uniquely determines the polynomial H,(u). If f(z) is an 
arbitrary polynomial in 2 then by (3.3) 


flut1+H) -Aflu+H) = Q-A)flu), 
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from which it follows that for given f(z) the difference equation 
(3.4) glut1) -Ag(u) = (1-a)flu) 
has the unique solution 
glu) = glu|A) = flu+H). 
It follows at once from (3.2) that 


(3.5) H?(u) = nH,_ 


and generally 


= Ri (u) (n>), 


n-r)! 
which implies 


n 
(3.6) H(utv) = (). 
r=0 


If we differentiate (3.2) with respect to A we get 

1 } 
(3.7) Hy + | A) = (u- 
which reduces to (2.12) for uv = 0. We also note that (3.2) implies 
(3.8) A") = (-1)"H,,(u, A). 


We remark that by (2.10) and (3.1) we have 


nr 
r=0 s=0 


Again, if we put 
(3.10) tty m) = 1-9 
j=0 
then it follows from (3.3) that 
1-A 


The polynomial /,(0, m) is usually called a Mirimonoff polynomial; a more 


(3.11) f,(u, m) = 
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general polynomial is discussed by Vandiver [14], see also Bachmann 
p. 117]. 
To get a multiplication theorem for H,(u) we consider 


m—| 


m—i—r ( rj _ m—1—r 


r=0 


which yields 


(3.12) m” | a”) = LATH (mu). 
m 1-A 
f= 


An interesting special case of (3.12) is obtained by taking A = ¢, where 
Then (3.12) reduces to 


m—-1 
(3.13) a” CTH, (wf | é) =H (ul 
r=0 
Nielsen [9, p. 54] has proved that the multiplication theorems for the 
Bernoulli and Euler polynomials characterize those polynomials. Suppose 
now that f(u|A) is a polynomial in wu of degree n that satisfies the equa- 
tion 


m—1 
ri \m 
r=0 


for some value of m>1. Put 


n 
= 
s=0 


Then by (3.12) and (3.14) 


= 
1-A ni 
n=0 


MATHEMATICS MAGAZINE (May-June 


n n 
Aga" (mu | a) = A LAH (mu |. 
s=0 s=0 


This requires that 
(3.15) [A™] =A,[A] (O<s<n). 


Now assume that (3.14) holds for two values of m>1, say m, and m,. 
Then it is clear that (3.15) becomes 


A,[A]=0 (O<s<n). 
Therefore 
flu| A) = A | A). 


For (3.13) the situation is somewhat simpler. If g(u| 2) is a polyno- 
' mial in uv satisfying the equation 


m—1 
(3.16) m” > (uk | = g(mu| 2), 
r=0 


where 1, 1, then if we put 
n 
dul d=) (ul 
s=0 


and assume that (3.16) holds for one value of m>1, then it follows readily 
that 


where A_(¢) is arbitrary. 
It may be of interest to mention also an addition theorem satisfied by 


H,(u| A). Since 


2 
1-A 1+A , _ _1-A 


e7-v “4A 


it follows at once that 


n 


r=0 


| 
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We note also that from the identity 


(1-A)+A) _ _ 2AC-A?) 


follows 


(3.18) (1+A)H(w| (uw | = 277 (4 a‘) 


while from 


follows 


n 
(3.19) (2A-1) (M)H AH, _ (| = AH, | A) -(1-A)H (ur | 1-A). 
r=0 


4, It is familiar that the Bernoulli polynomial B (uw) may be defined 
by [10, Chapter 2] 


(4.1) 


’ 


where ¢ is a primitive m-th root of unity and m+r. This evidently implies 


s=0 ¢ 


For m|r, on the other hand, we have the multiplication theorem for B Ma): 


(4.3) 


255 
x” 
n=0 
then 
oo 
n ru 
n=0 s=0 s=0 
net" 
n= 8 
m B,, (+i) = B, (mu) 
s=0 
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Multiplying both sides of (4.2) by ¢%, summing over r and using (4.3), we 


get 
m-1 


(4.4) ui = B, (mu) +n 


r=] 
where 0<t<m. 
We recall that the Bernoulli function B,,(u) is defined by 
B(u)=B,(u) (O<u<i), B(u+l) = B 
Similarly we define 
(4.5) Hy (uel | = CH (ul 2), 


where ¢ is some root of unity # 1. With these definitions of Hu |) and 
B (2) it is easily verified that the formulas (3.8), (3.12), (3.13), (4.2), 
(4.3), (4.4) hold for the barred functions; in particular in (4.4) the restric- 
tion 0<¢<m is no longer necessary. 


We remark that for m even and 2" = ~1, (4.2) reduces to the known 
formula 


m” (-1)°B,, (+4) =- (mu), 
s=0 


where E(u) is the Euler polynomial of degree n-1. 
For €" = w, where w?+ +1 = 0, m = 3, u = 0, (4.2) becomes 


For n even, it is known that [10, p. 22] 

B,(4) = B,(2) = 
from which it follows that 
(4.6) = (3"-1)B, (n even) . 
On the other hand for n odd > 1 we get 
(4.7) nH, lo] = 3"0B,(4) 


Again for ¢" = i, m = 4, we get 
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It follows that 
(4.8) ni-1)H, _ = 212"-1)B, (n even), 
(4.9) ni+ DH, _ = 2B (4) (n odd > 1). 
For \ = —1 we.have [10, p. 28] 
B. 
(4.10) H,,_,[-1] = = 201-2") =. 


5. We now obtain some congruences satisfied by H,,. If in (2.8) we 


take f(z) = 2™(1-z2”)", we get 


n 


s=0 
Assume that p is a prime such that 
(5.1) (p-1)p°—* | w; 
then by Fermat’s theorem 
= 0 ( 


We thus obtain (Frobenius) 


k 
s=0 


valid for k>0, n>0, provided w satisfies (5.1). This result is referred to 
as Kummer’s congruence for H,. In (5.2) A may be an indeterminate or an 
algebraic number such that (p, 1-A) = (1). In particular A may be an /-th 
root of unity, where / # pf ; 


By means of (5.2) it is proved in [5] that the coefficient A, , of (1.3) 
satisfies the congruence 


(5.3) = s (mod p°*), 


where p/-!<s<pi, n>e, b = pit8-\p-1). 
Another interesting result, also due to Frobenius, is 
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(mod p°), 


(5.4) 


where again w satisfies (5.1). Indeed we can prove a slightly more gen- 


eral result by using (3.11) with n = w, u = 7, m = 0 (mod p°). We get 


m-1 m-1 


=O 


A-1 
where the integer r, is defined by 
(5.5) r=-r, (mod p) (0<r,<p). 
We have finally 


_ 1_)P-1-ro e 
(5.6) H fr) = 1-A (mod p°), 


which reduces to (5.4) for r = 0 
The interesting congruence 
(5.7) = _ (mod p) 


is due to Mirimanoff [8]. It can be proved rapidly as follows. From (3.11) 
we get 


2 
(1-A)R ,_ = (j41) 
j=0 


so that 


d 
5.8 ofA 


Replacing A by 1-A, we get also 


R, 
= 


Consequently 


p-2 


THE TREE 
OF 
MATHEMATICS 


IMPORTANT FEATURES 


Presents epitomes of twenty-five of the main branches of 
mathematics, 

Starts on ground familiar to almost everyone. 

Is written in simple language except where the very nature of 
a topic requires specialized terminology. 


Problems designed to broaden and fix the concepts are supplied 
in the first seven chapters—algebra through calculus. 


- Puzzlements are omitted in the interest of brief presentations of 
basic principles. 


Extensive bibliographies are given for the newer courses such as the 
theory of games and dynamic programming. 


IMPORTANT USES 


Home study for those not taking regular class work. 
Broadening one’s mathematical horizon in a minimum of time. 


Reference book for anyone working in mathematics or its 
applications, 


A text book for survey (general) college courses in mathematics. 
THE FIRST SEVEN CHAPTERS 
A cure fo: those who have the mistaken notion that they “just 


cannot get mathematics”. 


Preparatory reading for students planning to specialize in mathe- 
matics, physics, engineering or other courses using mathematics. 


A text book for general courses in high school (might be time 
for additional exercises). 


Library of Congress number 58-847) 


From the Preface to: 


“THE TREE OF MATHEMATICS” 


There is a great deal being published about mathematics these 
days, and that is fine; but this book is mathematics in the sense 
that it presents the epitomes of the main branches of the subject 
beginning with high school algebra and extending far into graduate 
work. 

When writing this treatise the authors gave great attention to 
making it both meaningful and understandable, an art in which most 
of them are pastmasters. The practice was to start on ground familiar 
to everyone and construct a highway, free from road blocks, through 
the wonderful world of mathematics. 

This book is, for the most part, a response to requests for 
source material from two classes of people: those who need an ever 
increasing knowledge of mathematics in their jobs, such as en- 
gineers; and others who have gone little if any beyond arithmetic 
and either need more mathematics or just want to know “what it is 
ali about.” Inherent in a satisfactory response to these requests 
is an answer to the needs of teachers and students of mathematics 
and related subjects who desire to extend their horizons by home- 
study, and to do so in a minimum of time. 

In terms of class-work, the first seven chapters and perhaps 
selections from later ones would constitute about a one semester 
survey course in high school, while by passing swiftly over the 
first three or four chapters the entire volume could be covered 
in a two semester survey course in college. 

But classwork is always slower than effective individual study. 
The latter could reduce the above periods to months or possibly 
weeks by study at home during the evenings. 

Due to its broad coverage, this treatise offers an excellent 
opportunity for “reading up” on specific topics. By looking a topic 
up in the index you will find references to the page or pages where 
it is discussed or its meaning depicted by usuage. Of course,a 
similar procedure is possible with textbooks, but it would require 
some two dozen of them to cover all the topics in this treatise. 

Gratitude is due the many authors who have taken time from 
their crowded days to work on this unusual project, especially 
Professor Edwin F. Beckenbach for his exceptional support and 
the name “The Tree of Mathematics,” and D. H. Hyers for his val- 
uable advice on editing the book. 


Glenn James, Managing editor. 
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where C is independent of A. Clearly C = 0 and (5.7) follows at once. 
Again by (3.11) 


p-2 
j=0 


so that 


[A] . 


d ‘ 


Then by (5.7) 
now using (5.8) and (5.9), we get 


so that 
(5.10) -3R = + (1-2)? (mod 7), 
where of course p>2. This result also is due to Mirimanoff. 
We note that from (2.1) and 
(?=1) = (-1)" (mod p) 
it follows that 


p-l 
(5.11) (-1)'H, = AH,_, (mod p), 
r=0 


p-l 
1 
(5.12) p> (A-1)H, -1 (mod 
r=] 


Finally from (3.11) follows 


p-l 
(5.13) (mod p), 
f= 


where p = 2k+1 and (r/p) is the Legendre symbol. 
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PROPORTIONAL METRICS IN V VARIABLES 
H. Randolph Pyle 


Introduction. 

Suppose that we are given a surface S with metric ds? = Edr?+ 
2F dedy+Gdy? and a surface S’ with metric ds’? = Ede’? +2F ‘de ‘dy’+ 
G‘dy”*. The point (2’,y’) of S’-is related to the point (2, y) of S by the 
transformation 2’ = $(z,y), y’ = W(z, y). The necessary and sufficient 
condition for the conformal mapping of S on S’-is ds’? = A*ds*. We have 
shown in a previous paper (1) that this is equivalent to a set of linear 
relations between the partial derivatives of ¢ and w that reduce to the 
Cauchy-Riemann equations for isometric parameters. It is the purpose of 
this paper to find similar relations when the metrics are expressed in 
terms of n parameters instead of two. This will be accomplished by the 
use of matrix methods. 


Conditions for proportional metrics. 


non 
We shall express the first metric as ds? - > > a; ,de,dx, and the 
t=1 


second as ds“? = b dy js and impose the condition that ds” = 


The y’s are related to the z’s by the transformation ome. 
a,), ( = 1,++n). We shall make the usual assumptions concerning the 


existence and the continuity of the partial derivatives of the ¢,, and that 


0d. 

the Jacobian does not vanish. Let the symbol ¢;, stand for eed so that 
k 


_* . The Jacobian J is the determinant of the matrix 


k k 
Ono 
Let 
261 
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(a; 5 = Ory bi; = b 5) be the matrices of the coefficients of the differentials 
inds* and ds “* respectively. Designate the determinants of these matrices 


by |A|, |B|, and assume that |A| 4 0, |B| 4 0. Let dX = (dz ,,dx,,++-dz,), 
dY = (dy,,dy,,+ dy,) be row vectors and (dX)", (dY)™ be the corres- 


ponding column vectors. 
Using matrix multiplication, we have ds* = (dX)A(dX)" and ds” = 


(dY)B(dY)", Since dy, = 2 (4Y)™ = and d¥ = (dX)AT. 
=1 


Then ds’? = (dX)A7BA(dX)7. When ds”? = A2ds?, ATBA = 
Let ij represent the cofactor of $3; in J, and 


in 


2n 


0 
0 


n 
where / is the identity matrix. Now let Li; = Pe b. 

=1 
and 


Since ATBA = 424, A(ATBA) = A(A2A), or (AAT)(BA) = A2(AA), and 
JL = Also |AT||B||A| = |A?4|. But |AT| =|A|=J, and 


2 
so that J?|B| = and A? = 
J A 1/n p(n—2)/n 
Since JL = \7L, if we call p= , we have L = pL, 


and = or 


n 
1 


k=1 k= 


Ono Pnn Pin Pon Onn 
Then 
«as 
AAT = | 
a 
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as the necessary conditions for ds’? = A7ds?. 
Since each of these steps is reversible when J, |A|, |B| # 0, these 


conditions are also sufficient. 


Special cases. 
1. When r= 


ds” = Ede? +2Fdz dz,+Gde? 


ds”? = + 2F ‘dy dy, + G 


EF 
_ [vy He 


gives us the relations previously derived by other 
methods (1). 


2. In an isometric coordinate system, aj; = b; 37 8 (Kronecker delta), so 


that |A| = |B) =1,n=J"%, = and 
When n = 2, these become the Cauchy-Riemann equations. 
Let W = (y,,-+-y,,) be a variable vector for which y; = ,, +++ ,). 


Then = W; = see The outer product of the n-1 vectors 


W W 
Po; ooo When 


j+v is a vector 2; whose components are 


But a; is orthogonal to the vector space spanned by the n-1 vectors men- 
tioned above, so that Ws is orthogonal to this space. Hence the set of 
vectors W,W,,-+W, is an orthogonal set. 
1. H. Randolph Pyle, Conformal Mapping of Surfaces. Duke Mathematical Journal. Vol. II, 


Whittier College, Whittier, California 


CARRY-OVER 


Colonel John W. Ault 


It is probably safe to say that every mathematics teacher believes 
that any student will derive great benefit from a broad background in 
mathematics even though his life work may be in a field which requires 
only a basic skill in arithmetic. This feeling is based on the fact that 
there is a certain amount of “carry-over” of the mental discipline to which 
he has been exposed. 

Without any attempt to open the whole subject of carry-over, since it 
is one that most people consider debatable, the following remarks are de- 
voted to the simple approach to problem solving. By the simple approach, 
we mean the basic procedure of isolating the question to be answered and 
listing the given facts which can be brought to bear upon the question. It 
is far too evident that most students believe that (a) Given, (b) To Find, 
and (c) Solution approach to a problem should be left behind at the com- 
pletion of high school plane geometry. 

It is believed that if the student can be made to idealize this basic 
approach to problems, whether they be in mathematics, physics, econom- 
ics or whatever, that there will be some carry-over. In a military school 
it is normal to require that a student’s room be kept in apple-pie order at 
all times, but if he does not build up in his mind the fact that neatness is 
a better way of life and idealize this trait to the point where he is slight- 
ly dismayed by others who do not share this belief, then his first bachelor 
quarters after graduation are pretty apt to be sadly disorganized. 

Just how the teacher can build up this idea to the ideal level is 
hard to say. Of course, great benefit can be derived from constant em- 
phasis of the principle in lectures and drill sessions, and grasping every 
opportunity to show the student who has difficulty with a problem that he 
would have avoided his difficulty by using this approach. Perhaps in the 
long run the most effective teaching device will be the fact that the 
teacher himself applies this method to all problems met in mathematics 
or elsewhere. 


United States Air Force Academy 
Colorado 


A NOTE ON SIMPLE CORRELATION 


Maurice Fréchet 


Translation from the French 
by 
C. de la Menardiere 
Department of Modern Languages 
University of Alabama 


In Volume 22, No. 2, 1958, pp. 57-69 of Mathematics Magazine, ap- 
peared an article by C.D.Smith entitled, “On the Mathematics of Simple 
Correlation.” The aim, as given in the article, is to state the problem as 
it usually occurs and show how the measures of correlation may be de- 
rived by elementary algebraic calculations. This aim was completely ful- 
filled. The purpose of this Note is to add some comment regarding inter- 
pretation of measures of correlation which follow from theoretical con- 
siderations. 

The paper states, p. 58, “The method of correlation is designed to 
measure the strength of common influences between two variables.” We 
may note first that the coefficient of correlation r does not always indi- 
cate what influences may exist. Yule, and others, (1), have found cases 
where the measure has a value near its maximum where no common in- 
fluence is known to exist between the two variables under consideration. 

In fact, the thing that can best measure the coefficient of correlation 
is the thing that can measure a good index of correlation. We seek the 
total accuracy of the determination of one of the two variables when the 
other is given. In other words, an index of correlation has nothing to say 
as to the causes of this accuracy or the lack of accuracy, not even the 
existence of the causes. It gives a numerical fact. For instance X, is 
more or less determined by knowing X,. This is why the term, “nonsense 
correlation,” seems deceiving to us. fe all cases, the value of an index 
of correlation has a meaning if the numerical data from which it has been 
extracted have been correctly taken. But that meaning is only one of the 
characteristics of the numerical table of correlation. 

When a good index is equal to 1, we may state four conditions, name- 


ly: 


I. To all values of X, there corresponds only one value of X ,. 
II. Conversely when it is equal to 0. 
III. X , is constant when X, varies. 
IV. And conversely. 
Moreover it is known that the classic coefficient, r, verifies only two of 
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these conditions (I and IV) (1). In some cases it is a bad index, the use 
of which must be controlled. The value of a good index of correlation 
does not dictate to the statistician the above conclusion. It invites him 
to look for the causes of the mathematical dependence, when for instance, 
the index is close to one. 

Consider the coefficient of linear correlation given by the mathe- 
matical formula: 


2 
No 


(1) r 


Also consider the correlation ratio, 7 as defined by n = —. 


1 
(2) We ask: 1. Why is it a correlation measure? 


2. Why 7 is introduced if r is already known? 
Karl Pearson, having noticed the faults above concerning r, tried to 
make a new formula. 
A calculation, quite elementary, permits one to see better the pass- 
ing from r to 7 (1). 


Indeed if one has p = £, one has according to (1) and (2) 


n ? 
p = . 
Then since 22, = 0, 


No, 4 


where m is the average of X, for 2, = x, and te is the variance of m. 
From this it follows that: 
r= pn 

Here p is simply what becomes r when one replaces the set of points 
X,, X, by the line of average X ,, when X, is given. Then p is, like |r|, 
at the most equal to 1 in absolute value and one has |r| < 7 < 1. One sees 
that p does not vary when the set of points (X ,, X,) stretches or expands 
around the line of average. \ 

Consequently the coefficient r contains a strange factor, p, Which 
tells us nothing about the proposed sum. Is the set spread or contracted? 
Are the values of X, more or less determined when one knows \2? By 


eliminating that strange factor, one has this time an index 7, which veri- 
fies three of the conditions (I, II, III, IV). And even if it does not verify 
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exactly the condition III, which is to say, if when it equals zero, X, is 
not necessarily constant, at least the average of X , is constant, for given 
X,, when X, varies. 

Sometimes one considers that 7 is much more sensitive than r to the 
grouping of data by intervals and some uncertainty arises from grouping. 
In my opinion, this fact indicates superiority of n over r. Because the 
grouping by intervals modifies not the real correlation between X , and X, 
but that which is represented on the correlation table. An intex of corre- 
lation should reflect this change and consequently warn the statistician 
not to much reduce nor to much spread the grouping intervals. 

It would be useful to add that, in spite of the superiority of 7 over 
r, neither r nor 7 will satisfy the four conditions above. Other authors, 
pursuing the improvement realized by Karl Pearson, have defined several 
correlation indexes satisfying the conditions I, II, III, IV. One such 

.measure for instance is the index called connection by Gini (2). It would 
be desirable, for practical purposes, that the conclusions drawn from the 

alculation of r and 7 by means of the calculation of the index of connec- 
tion of Gini be controlled. The numerical calculation was considerably 
simplified by Salvenini. 

We have explained in what manner r is inferior to 7 (itself insuffi- 


cient) establishing by an elementary and rapid calculation the formula 


r= pn. 
Here p is the linear coefficient of correlation that one could obtain re- 
placing all the values of Y, when X is given, by their average. It is a 
geometric characteristic of the curve of the averages, equal to 1 when 
that curve is a straight line, and, it is only in the calculation of 7 that 
the dispersion of Y intervenes when X is given. 

Since the coefficient r gives satisfactory results when the couple 
(X ,X,) satisfies the Laplace-Bravais Law, called normal to two varia- 
bles, one expects to go back to that case, by transformation of the scales 
of X, and X, separately, in order that each of the new variables satis- 
fies the, “second Law of Laplace,” called normal Law. We have shown 
(3) that even after this transformation, the intended result has not been 
reached; the new law of the couple (X ,, X,) modified is not necessarily 
the binormal law. When the frequency laws, called marginal, of X, and 


X,, considered separately, are given, the linear coefficient of correlation 
r of any tables of correlation having those margins remains included be- 
tween two numbers r,, 7,, which are generally not equal either to -1 nor 
to 1, but simply included between -1 and +1 (4). These values of r,, r,, 


of contrary sign, correspond to two monotonic functional relations & ,, Rk, 
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between X , and X,. And reciprocally, if a monotonic functional relation 


exists between X , and X,, it is identical to #, or R, (5). 
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SOME FURTHER NOTES ON THE THEORY 
| OF CORRELATION 


C. D. Smith 


The paper, referred to in the Note of Fréchet, proceeds from the Re- 
gression Model as given by Yule, (1). The Note of Fréchet adds comments 
on the measures r and 7, listing certain limitations which follow from 
theoretical developments. In this note we add some comments which fol- 
low from the Regression Model, (1). Some additional references are given. 

1. It seems that Karl Pearson first used the Normal Probability Sur- 
face for two variables as a model, and defined the correlation r as the 
product moment of deviations of the variables from their respective means, 
with the standard deviations as units of measure. The two straight lines 


located on the plane of (X ,, X,,) by the method of least squares were used 
by Yule as the corresponding regression model for representing r, (1), 


page 188. The product formula for the frequency model is r = ———. Yule 
1%9 
from the regression equation z, = 


le 


gives the equivalent formula r = 


12, 


2. When the points in the z,-arrays spread widely about the corres- 
ponding points of the regression line, the value of r is relatively small, 
and estimates of 2, are not materially improved by using linear regression. 
When points in the arrays cluster closely about the corresponding points 
of the regression line, the value of r is relatively large, and estimates 


of x, may be improved by using linear regression. 


where o,. is the standard deviation of the linear estimates. 


le 


3. The correlation ratio 7 is not less than rfor a set of values (X ,, X ,) 
because the value of 7 is determined by the manner in which values within 
an array vary about the mean of the array, and the dispersion about the 
mean is less than that for any other point of an array. The values of r 
and 7 are equal for a given set (X ,, X,) if the means of arrays fall on the 
regression line. The values of r and 7 are abstract numbers calculated 
without regard to the manner in which values of X, and X, may be paired. 
Consequently they give no information regarding possible causes of cor- 
respondence. The worker must plan his problem so that some reasonable 
meaning may be assigned to the results. 
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4. Values of correlation may be relatively large in some cases where 
there is no apparent cause for correspondence. H. L. Rietz illustrated the 
extreme case where 7 is greater than r by using the function y = sing, 
0<a2<2n.Here r=0 and 7» = 1. The comparison illustrates a case of non- 
linear regression. 

5. Yule, (1), page 195, refers to correlation as a measure of the 
strength of the influence of one variable upon another. The measure does 
not identify the influence. 

6. R.A. Fisher, (2), page 190, states the case as follows: “The cor- 
relation between A and B measures, on a conventional scale, the impor- 
tance of the factors which (on a balance of like and unlike action) act 
alike in both A and B, as against the remaining factors which affect A 
and B independently.” 

The case for common factors may be illustrated by the following ex- 
periment. Toss 5 coins on the table and count the number of heads (X ,). 
Leave the 3 coins nearest the edge of the table and toss the other 2 coins 
down with them. Count the number of heads (X,). Repeat the experiment 
to obtain a sample of pairs (X ,, A.,). The two series of counts are corre- 
lated and we say the 3 coins common to the two counts of heads consti- 
tute the common factors. In performing the experiment there could be other 
factors such as imperfect coins, or bias in the manner of tossing them. 

7. Rietz, (3) pp. 80-81, gives a probability model for correlation in 
the form (2, y)dady = g(r)A(x, y)drdy. When A(z, y) is a function of both 
r and y, the two variables are said to be correlated. 


LIST OF REFERENCES 
(1) G. Undy Yule; Introduction to the Theory of Statistics, 6th ed. 1922. 


(2) R. A. Fisher; Statistical Methods for Research Workers. 


(3) H. L. Rietz; Mathematical Statistics, Carus Mathematical Monographs, 
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MISCELLANEOUS NOTES 
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Charles K. Robbins 


Articles intended for this department should be sent to Charles K. Robbins, De- 
partment of Mathematics, Purdue University, Lafayette, Ind. 


A FURTHER GENERALIZATION OF NEUSTADT’S LAW 


George Bergman 


In his article “A Mathematical Munchausen” in the Nov.-Dec. °56 
issue of this publication, Paul Schillo presented, amidst a large quantity 
of satire, a sequence of numbers called “Neustadt’s Base” and a certain 
conjecture, “Neustadt’s Law”, concerning this sequence. The sequence 
was: .2, -l, 2, —2, 1, —.2; and the “law” that if this series ,be multi- 
plied, term by term by any succession of arithmetic progressions, the re- 
sulting sequence will have zero sum. E. g.: 

2 2 -2 -.2 
x2 
-2 
x3 
-6 6 8 -15 
x1.80x1.44 x1l.08 x.72 x.36 


1.80-8.64 6.48 5.76 -5.40 


In the Jan.-Feb. °58 issue, W. W. Funkenbusch stated the correction 
that the number of successive multiplications must be limited to 4, and 
made some observations concerning multiplication by progressions of 
higher degree. 

However, all this is merely one case of a much more general law, 
with numerous applications, namely: 

If P is a polynomial of degree less than n, 


n 
=o 
i=0 
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The proof of this consists of three parts: 
a) A polynomial P of degree less than n can be written in the form P(z) - 


a(,,* ,)+@(z) where Q is a polynomial of degree less than n-1 
To show this, let A equal the coefficient of 2”~' in P(x) (which may 


be equal to zero). Now (,,*,) is a polynomial inz of degree n—-1, in which 


the coefficient of is . Thus, if we subtract (n-1)!A( ,) from 
n—1)! 


P(x), we will get a polynomial of degree less than n-1. I.e., letting 
(n—1)!A = a, we have 


P(x) ,) = Q(z) 


or 
P(x) = 7 ,)+@(z) of degree n-1) 
QED. 
b) | = 0 
2=0 


The only values of ¢ for which Pes # 0 within the indicated range 
are n and n-1l. 


Thus 


c) If for all polynomials Q(z) of degree less than n-1, 


n—1 
> = 0, 
7=0 


then for all polynomials P(x) of degree less than n, 


n 
Seni ”)P(i) = 0. 
i=0 

To show this we begin by applying a, and then b to the expression given: 


n n 
i=0 
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n n n 

7-0 i=0 i=0 
Now, as Q is a polynomial of degree less than n-1, we can prove this 


last expression equal to zero if we can break it into expressions of the 
form given in our hypothesis. This can be done as follows: 


n n 
n n 
t= 


The final term of the first of these sums is zero, and first term of the 
second sum is zero. We can adjust the ranges accordingly. To make the 
second conform with the form in the hypothesis, we shall then substitute 
j = t-1. Q(j+1) can be expressed as a simple polynomial in j, of the same 
degree as Q, which we shall call Q’: 


n 
t= 


n-1 n 
i=0 


= 0+ 
j=0 
j-0 


This proves c), which, together with the self-evident case when n = 0, 
constitutes an inductive proof of the theorem. 

Now “Neustadt’s Base” merely consists of the numbers (-1)( 5), 
disguised by division by 5. Since the product of less than five arithmetic 
progressions is a progression whose general term is a polynomial of de- 
gree less than five, the sum of the products of corresponding terms of 
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less than five arithmetic progressions, multiplied by the corresponding 
terms of Neustadt’s Base, will equal zero divided by five, or zero. (“I. 
Neustadt” did not believe in mathematical induction — butone can’t please 
everybody! The particular case of “Neustadt’s Base”, or any particular 
case, can be proved deductively by multiplying the numbers of this “base” 
successively by a, a+d, a+2d, ---; then by a4, a’+d’, then 

The applications, as I mentioned, are numerous. For instance, take 
the method commonly used for finding the degree of a progression {S_}. 
We find first the differences between successive terms, AS, and then the 
differences of the differences, A?S, and so forth, until we reach an n such 
that all the A*S equal zero: whereupon we can state that n-1 is the de- 
gree of the progression. This can be explained in terms of the above 
theorem when one considers the fact that the éth A’S is 


n 
> - 


Another use for it is in the summing of progressions whose general term 
is Since the numbers (-1)) are the coefficients of the in 
the expansion of (1-r)", if one multiplies the series by (1-r)¢+! where d 


is the degree of P(), then the coefficient of r” becomes 


d+1 

(-1)( 4+ 

t=0 
or 0, except for the first and last d terms, where the sum will be incom- 
plete for want of non-included terms. Thus, all but 2d terms will drop out, 
and dividing through by (1-r)*, one will have a closed expression for the 


sum. 
Example : Evaluate 


S = 1447+ Or? + (n—2) (n—1) 299-2 
Solution: Multiply through by 1-3r+3r?-r° to obtain 


(1-r)* 


*as, for instance, Proposal 333. r= = 


Stuyvesant High School, New York, N. Y. 


A NEW LOOK AT AN OLD PROBLEM 


Charles E. Wingo Jr. 


I have always been deeply interested in mathematics, and for the 
last ten years mathematics has been a source of recreation and enter- 
tainment for me. It was thru this that the old problem 


a*+y = ll, y?+2 = 7 


was revived. It was over fifty years ago that I first heard of it. What has 
stimulated interest in this problem is the fact that it has proved to be as 
much a puzzleas a problem. 

Equations of the form 


z*+y=a, 


were pronounced insoluble by Prof. Quimby of Dartmouth, except by Des- 


cartes formula for solving bi-quadratics. 
Those who are further interested in some of the solutions to the 


problem 
y?+2=7 


are refered to the Am. Mathematical Monthly, Vol. V, pages 291-2-3-4-5 
and Vol. VI, pages 13-37, where will be found about eight solutions, 
some of which are quite complicated; bi-quadratics, factoring, and other 


methods. 
In the problem as treated here, an effort is made to develope a method 
of solution whereby all problems of the form 
y?+r=6 


can be easily solved, where z > y and az, y, a and 6 are positive inte- 
gers, provided of course that integral solutions exist. This means that 
arbitrary values of a and 5 can not be taken, as 
The method derived here is general in scope, subject to the limi- 
tation mentioned above. Consider the two equations 
and 
y?+2=b. 
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We have 

=a-b 
or 
(1) (x-y)(2+y-1) = a-b 


In my study of this problem over the past few years, my interest centered 
on equation (1) above. For the solution, if any was to be found, would 
hinge on this equation. For it is obvious that (a—5) is the product of two 
factors (x—y) and (x+y -1), and the key to the solution was to find which 
of these equate to (x-y) and which to (r+ y-1). 

At this stage the problem rested for some time. Then some time later 


it flashed to me that Jarry and Vytte could be written 
= Va 
= b= ytd 
Va-yb = #-y+ce-d 


Since z > y, we know that in taking the square roots of x7+y and y?+z2, 
the amount that c adds to z is less than the amount that d adds to y. 


Therefore d is larger than c. So a ~ Vb will be less than the true value 
of z-y by the amount of d-c. So the value of the factor of a—4 to equate 


to z-y will be the integral value next greater than Va-vb and (r+y-1) 


would be equal to the other factor. 
To show how the method is applied to problems we consider the 


following. 


Factors of 4= 1-2-2 
«aa 
= S87 
Va-yb = .65 
The first integral factor larger than .65 is 1. So 
r-y=1 
(continued on page 288) 


PROBLEMS AND QUESTIONS 
Edited by 
Robert E. Horton 


Readers of this department are invited to submit for solution problems believed to be 
new and subject matter questions that may arise in study, in research, or in extra-aca- 
demic situations. Proposals should be accompanied by solutions, when available, and 
by any information that will assist the editor. Ordinarily, problems in well-known text- 


books should not be submitted. 
Solutions should be submitted on separate, signed sheets. Figures should be drawn 


in India ink and twice the size desired for reproduction. 
Send all communications for this department to Robert FE. Horton, Los Angeles City 


College, 855 North Vermont Ave., Los Angeles 29, California 
PROPOSALS 


376. Proposed by C.W.Trigg, Los Angeles City College. 

Identify the following unorthodox definitions: (1) the earthy or stony 
substance in which an ore or other mineral is bedded, (2) fully sufficient, 
(3) the four brightest components of @ Orionis, (4) type of speech fanci- 
fully exaggerated, (5) oddity, (6) a Franciscan friar, (7) a deposit formed 
in a liquid vegetable extract, (8) a large molding of convex profile, (9) 
fancied, (10) the series of air bubbles made by the breath from an otter 
under water, (11) a drain to carry off filthy water, (12) one sixteenth of a 
fluid drachm, (13) scheme, (14) depicted as walking, (15) absent-minded, 
(16) seekers, (17) complete, (18) invalid, (19) an interpreter of music. 


377. Proposed by J.M.Gandhi, Thapar Polytechnic, Patiala, India. 
Prove that 
3n+] 
3n+1 an+2 
7=0 
when n is of the form 2m, 4m+3 or 4m—1 respectively. 
378. Proposed by Barney Bissinger, Lebanon Valley College, Pennsyl- 


vanta. 
Let a,, ;a, be arbitrary and n(ntl)a,,, = n(n-l)a,-(n-2)a,_, for 


a>0. Pd Se... 
0 nm 


379. Proposed by J.M.Howell, Los Angeles City College. 
Two men fight a duel. They both fire at a given signal. If both are 
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alive, they fire again at a given signal, and repeat the process until at 
least one of them is dead. If the probability that A kills B on any round 
is p, and the probability that B kills A on any round is +, find the proba- 
bility that A is alive, B is alive or neither is alive after: 1) n rounds, 
2) after an infinite number of rounds. 


380. Proposed by Huseyin Demir, Kandilli, Eregli, Kdz., Turkey. 
Solve the system of equations 


a(2—a) + u(a+u) = 0 
yla—b) +uly+u) = 0 
—aly-c) +ulz+u) = 0 


where abc #0 and = y7 1. 
381. Proposed by George M. Bergman and Melvin Hochster, Stuyvesant 
High School, New York. 

Show that if the circumsphere and the insphere of a tetrahedron are 
concentric, four spheres can be drawn each tangent to each edge of the 
tetrahedron, extended if necessary. 


382. Proposed by C.N.Mills, Sioux Falls College, South Dakota. 


2/3 


For the hypocycloid 2?/*+y?/% = a?/3, («, 8) is the center of curva- 


ture. Show that «+f = (a 
SOLUTIONS 


Late Solutions 


350. J.L. Brown, Jr., Ordnance Research Laboratory, Pennsylvania State 
University; M.Morduchow, Polytechnic Institute of Brooklyn. 


The Wine Merchant 


355. [November 1958] Proposed by P.B.Jordain, New York, New York. 

A wine merchant had a small cask containing 100 gallons of wine. 
In order to make more money, he decided to replace each gallon he took 
out of the cask by a gallon of water. This he did n times. Finding that 
he was losing customers, he naively tried to undo what he had done by 
selling his watered wine from the same cask, but replacing each gallon 
sold of the mixture by a gallon of pure wine, feeling that by this method, 
at the end of an additional n gallons of watered wine sold, he would have 
a cask of full wine again. He was unfortunate enough in having hit upon 
a total number of gallons sold such that the wine contained in the cask 
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was at a minimum at the end of this operation. The question is, how many 
gallons of watered wine had he sold in all? Assume that the merchant 
dealt only in full gallon sales. 

Solution by David F. Finlayson, Worcester Polytechnic Institute, 
Massachusetts. After n sales the fraction of wine remaining would be 


given by (a. The gallons of wine remaining would be given by 
100( =) = am The gallons of water left would be given by 100- 


ane After repeating this process n more times by replacing the mix- 


ture sold by pure wine instead of water as in the first n sales, the gallons 


of water left would be given by (98 "(100 - GOI. The wine in the 


cask would then be given by 


_( 99 yn 
100 [100 G. 


or 
an 
G = 100-99)" _ , _(99) 
(100)"~* 
G = 100 -(99)"(100) + (99)?"%(100) 
Differentiating with respect to n and setting a = 0 gives 
n 


=—log,99+ log, 100+ (99)"log ,99(100)~(2) + (100)~log ,100(99)"(-2) = 0 
n 


log,99- log, 100 = 2( log,100) 
1 = 2(-29)" 
100 


0.5 = (0.99)” 
n = 69 * 
2n=138 2n-1=137 


Since the first gallon sold was pure wine, the wine merchant sold 
2n-1 = 137-gallons of watered wine. 

Substituting the value of n in the expression for the number of gal- 
lons of wine left, we find that the wine merchant ends with 83.57 gallons 
of wine in the cask. | 
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*It can be immediately seen that (.99)®® can not equal exactly 0.050 
since the last digit of .99 to any power can not be 0. Solving the equa- 
tion by means of a slide rule or logarithms produces an answer very close 
to 69. It is better, however, to square both sides of the equation and 
solve for 2n. By doing this an irrational answer between two succesive 
integers will be found and of course the even integer is the answer de- 
sired. The value for 2n is then found to be 138 and n must equal 69 and 
not 68 or 70. 

Also solved by Stephen A. Andrea, Oberlin, Ohio; Philip Fung, Idaho 
State College; Joseph D. E. Konhauser, Haller Raymond and Brown, Inc., 
State College, Pennsylvania; Sam Kravitz, East Cleveland, Ohio; Richard 
Mittleman, Los Angeles City College; Lawrence A. Ringenberg, Eastern 
Illinois University; C.W. Trigg, Los Angeles City College; DaleWoods, 
Idaho State College; and the proposer. 

Several solutions were received with the value of n incorrectly ap- 
proximated by 68 instead of 69. Ringenberg pointed out that this problem 
is a special case of problem E962 appearing in the March- 1952 issue of 
the AMERICAN MATHEMATICAL MONTHLY. 


A Reciprocal Sum 


357. [November 1958] Proposed by Joseph Andrushkiw, Seton Hall Uni- 
versity, New Jersey. 
Show that if d>0, 


1/(k?-2dk+2d?) = k+30/4d, O<k<1/2d? 
k=0 


Solution by Arne Pleijel, Trolhattan, Sweden. Let the proposal be 
restated in the form: 


2_Odk+2d2 4d = 2d? 
k=0 


Then we have 


dx 
Jo 


dx 
24d? 


30 oo 
4d 


PROBLEMS AND QUESTIONS 


k+1 


>» maha 


— Sa is decreasing in z, we have 
(x-d) +d 2 


and thus, as 


oo 


1 


and so the proposition is proved. 
Also solved by the proposer. 


The Bouncing Ball 


358. [November 1958] Proposed by C. W. Trigg, Los Angeles City Col- 
lege. 

A ball having fallen from rest a vertical distance A, strikes a stone 
protruding from a wall and bounces off horizontally without spinning. If 
the distance of the stone from the ground is s and the coefficient of res- 
titution is e, show that 

a) the ball will strike the ground at a distance 2\/she from the foot of 
the wall; 

b) the inclination of the stone to the horizontal is arctanye. 

Solution by Philip Fung, Idaho State College. Let u, v be the veloci- 
ties of ball before and after impact and their components parallel and per- 
pendicular to the surface of stone be u,, u, and v,, v, respectively. Let 
the inclination of stone to the horizontal be 6, the time taken by the ball 
to strike the ground after impact with stone ¢, and the horizontal distance 
covered by the ball during this time be PF. 

Now, the velocity of ball just before impact with the stone is u = 


vy 2gh. So 


u, = uSin 6; 


uCos @. 


Since u, is parallel to the stone surface and thus unaffected in magnitude 
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2d? (k-d)24d2.= 4d 

k=0 
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and direction by impact, therefore we have, 
v, = u, = uSin 
= eu, = eu Cos 6. 


But, the direction of v is horizontal. Thus, 
Tan @ = = 


From this we have, 


6 = arctanVe. 
Now, 
s = 1/2 gt? 
or, 
t= 
And, 


R=vt=t 


Hence, by substituting the values of ¢, v,, Vy and then 6 into the equa- 
tion above, we can thus solve that R = 2\/she. 


Also solved by William E. F. Apphun, St. John’s University, New 
York; Joseph M. C. Hamilton, Los Angeles City College; Joseph D. E. 
Konhauser, Haller, Raymond and Brown, Inc.j State College, Pennsyl- 
vania; Dale Woods, Idtho State College; and the proposer. 


(May-June 
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The String-of-Pearls Polynomial 


359. [November 1958] Proposed by Norman Anning, Alhambra, California. 

Prove that the string-of-pearls polynomial z* +1 can be expressed in 
at least one non-trivial way as the sum of two squares, if k is any even 
positive integer different from 2, 4, 8, 16+++. 

Solution by William M. Sanders, Mississippi Southern College. Let 
k be the n-th term of an arithmetic progression with first term 0; i.e., & = 
(n-1)d. Impose the condition that d be the maximum divisor of & such that 
n = 2p > 4 to obtain the progression 0, d, 2d, -+-, (n—-2)d, k. This is pos- 
sible since & is even and different from 2, 4, 8, 16 --». Now for k = 2m, 
d = 2q, and n = 2p, 


2p—2 
(-1)"4 min(v, 2p—v—1) 2M—20q 
v=1 
2p—2 
(-1)"4 min(v, 2p-v—1) 2 


v=l 


p-l 2 p-l 2 


v=1 v=] 


also svulved by L. Carlitz, Duke University; Arne Pleije!, Trolhattan, 
Sweden; and the proposer. 
A Difference Equation 


360. [November 1958] Proposed by Chih-yi Wang, University of Minnesota. 
Considering the higher differences of (”) with respect to z show that 


Solution by Arthur E. Danese, Union College, Schanectady, New 
York. In Rodrigues’ formula. for Jacobi polynomials [Szego, ORTHOG- 
ONAL POLYNOMIALS, New York, 1939, p. 66, 4.3.1]: 


2"! dx 


replace zx by ¢-1, n by n-2, B by z, « by r-n+m, and evaluate at ¢ = 1 


to obtain 
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which equals 


since 


(-1)°P 


The use of the explicit evaluation of the Jacobi polynomial : 
n 
=0 
{ibid, p. 67, 4.3.2] yields 


n—2 v+my m n 
v=0 


which is equivalent to A”("). 
Also solved by the proposer. 


A Variable Nine-Point Circle 


361. [November 1958] Proposed by N. A. Court, University of Oklahoma. 

A variable triangle inscribed in a rectangular hyperbola has a fixed 
vertex and the opposite side moves parallel to itself. Show that its vari- 
able nine-point circle passes through two fixed points. 

Solution by J. W. Clawson, Collegeville, Pennsylvania. It is well 
known that the orthocenter of a triangle inscribed in a rectangular hyper- 
bola lies on the hyperbola. Also that the center of the hyperbola lies on 
the nine-point circle. (C. Smith, Geometrical Conics, pages 165, 166.) 

Let A be the fixed vertex, BC the variable side which moves parallel 
to itself. Let the perpendicular from A to BC intersect the hyperbola at 
H. Then the fixed point H is the ortho-center of the variable triangle ABC. 
Since the nine-point circle of a triangle passes through the mid-point of 
AH, the variable nine-point circle passes through this point and also 
through the center of the hyperbola. 

Also solved by Huseyin Demir, Kandilli, Eregli, Ka., Turkey; 
Sister M. Stephanie, Georgian Court College, New Jersey; and the pro- 


poser. 


1959) PROBLEMS AND QUESTIONS 


Comment on Problem 337. 


337. [March 1958 and November 1958] Proposed by Victor Thebault, 
Tennie, Sarthe, France. 
Comment by William E. F. Appuhn, St. John’s University, New York. 
In order to eliminate having sides less than or equal to zero, condi- 
tions should have been placed as the solution was developed : 
1. In equation (1) 2 = 2pq, the condition pq > 0. 
2. In equation (2) y =.»?~ q?, the condition p? > p?. 
3. In the equations for p and qg we do have to place restrictions on 
the integers g and h. They are: 
(a) 940, (b) h £0, (c) g? # h, (d) gh > 0 if g? > A?, and (e) gh < 0 
if g?<h?. 
. In the final set of equations, we have to place the restriction 
gh>o. 
The above conditions will be satisfied by gh > 0 and g? > h?. 
>. Noting that, by equations (1) and (3), z and z are both even hence 
y is also even, a much quicker and easier complete solution may 


be obtained as follows: 
Let (1) @ = 4uv, 6 = 2!u?-v?| and hypotenuse e = 2(u?+v?) = 2d 


where u and » are different positive integers. In order that u?+v? = d? let 
(2) u - = |w?-y?|, 0. d = 27 +y7, where y are integers. Substi- 


tuting (2) in (1) to obtain the immediate solution: 
(3) a = Bry |x?-y?|, b = and c = Xx? +y7)?. 


QUICKIES 


From time to time this department will publish problems which may be solved by 
laborious methods, but which with the proper insight may be disposed of with dispatch. 
Readers are urged to submit their favorite problems of this type, together with the elegant 
solution and the source, if known, © 


Q 247. For what values of z is m?+n*-a?—b? > (mn-ab)x where 0 <a<m 


and [Submitted by M.S. Klamkin} 

Q248. Show that if & is any real number, the lines 24+hkr°y-62?y?- 
key? +y* = 0 cut 27+y? = 1 into eight equal parts. [Submitted by Norman 
Anning] 


Q 249. Show that for the smallest values of the angles, 2 arccos/1/3 = 
arc cos (1/3) = Il. [Submitted by C.W. Trigg] 


Q250. Prove that 19'° can not be represented as the sum of a fourth 
power and a cube. [Submitted by D. L. Silverman] 
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Q251. In a rectangular coordinate system (2, y), y = f(x). If the coordi- 
nates are transformed into new ones, (z‘,y’), by a rotation through a posi- 


dy’ d a? . 
tive angle «, express and in terms of and [Submitted by 
Henry E. Fettis] 
Q 252. Find the sum of 


[Submitted by M. S. Klamkin] 


Answers 


‘(udu ) = WINS 047 GARY OM 
= 


ay} JO Sepis UO Jo Selotjjooo ug 


Kp P 4 p q AT! — 


aM JO OY} JOPUN SI oy} BOUTS 


pus SIxe YIM g UB SeyvU (x)/ = A OY OUT] C4) J] 

= 619 = 6l IN| “(EL POW) 2 Idoeoxe 07 oq Avu ,A+ 
SNYL “(ET Pow) 6 40 = ,4 pow) ZT JO ‘T ‘0 = 
= [ = D Sasvo 

ul *2 = ¢/(2Z-9) sod +9 /D/\ Sod = [J SO9 


+ 9048 Og = LI) + UL 


PROBLEMS AND QUESTIONS 


> 10,09 < g OM 
Og SSE] OQ YQ Jopslo Ul SOD = Aq pus 
299 9S = pus = Gy ‘Y= OV “w = gy 
9 S09 <Q SOD UWZ— BABY OM “9 = 
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(A New Look at an Old Problem — Continued from page 276.) 


r+y-l=4 
22=6 2r=3 
2y=4 y=2 


Next take: 
2? +y=1052=a 
= 816=6 
a-b = 236 
Factors of 236 are 2.2.59 = 4x59. 
Va = V1052 = 32.43 
Vb = 816 = 28.55 


Va-yvb = 3.88 
The integral factor which is next and larger than 3.88 is 4. Then 
z-y=4 
z+y-1=59 


22=64 32 
2y = 56 y= 28 
For a further consideration of the problem, take the case where 


a= 6, Then: 
=a 
(x-y)(e+y-1) =0 
Since z and y are both positive integers, (r+ y—-1) is also a positive inte- 
ger. Therefore (r-y) is equal to zero. Therefore x = y. So the equation 
above becomes : 
a(z+1) = a 


It is seen from the equation above that a is the product of two consecu- 
tive integers. Hence a has such values as 1x2, 2x3, 3x4 etc. 


Amelia, Va. 
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